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1 .0 ABSTRACT 


This report contains the information needed to run a computer program 
for the calculation of the potential flow about arbitrary three-dimensional 
lifting configurations. The program contains a geometry package which 
greatly reduces the task of oreparing the input data. Starting from a very 
sparse set of coordinate data, the program automatically augments and redis- 
tributes the coordinates, calculates curves of intersection between components, 
and redistributes coordinates in the regions adjacent to the intersection 
curves in a suitable manner for use in the potential -flow calculations. A 
brief summary of the program capabilities and options is given, as well as 
detailed instructions for the data input, a suggested structure for the pro- 
gram overlay, and the output for two test cases. 
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3.0 INTRODUCTION 


This report describes the use of a computer program for the calculation 
of the potential flow about arbitrary three-dimensional lifting configura- 
tions. It is intended primarily as a user's manual for a geometry package 
which has been incorporated into the program, but input instructions for 
the complete program are included. The theoretical foundations of the po- 
tential flow method are described in reference 1. The geometry package, 
described more fully in reference 2, can significantly reduce the time and 
expense involved in using this potential -flow program. 

Input to the computer program consists of a set of flags to control the 
overall program execution, coordinate data to define the configuration, and 
flags to control the operations performed by the geometry packaqe. The 
coordinate data may be input using basically the original input format de- 
scribed in reference 3 (referred to as "Douglas inout") or the fo.'mat of 
the program described in reference 4 (referred to as "Deriv input"). The 
details of both these types of input format have been modified slightly to 
facilitate the incorporation of the geometry package into the program. 

The geometry package allows the defining coordinate data to consist of 
a yery sparse set of points, often with an order of magnitude fewer points 
than required for the actual potential flow calculations. Isolated compo- 
nents, such as wings, fuselages, etc., are automatically paneled, using one 
of several possible point distribution algorithms. Curves of intersection 
between the components are calculated and then components are repaneled so 
that adjacent elements on either side of intersection curves line up in a 
satisfactory manner for the potential flow calculations. 

Section 4.0 discusses, in general terms, the alternative capabilities 
included in the program. Further explanations can be found in references 1 
and 2. Section 5.0 describes the input data preparation in detail. Section 
6.0 suggests the structure of the program overlay to reduce the storage 
requirements to a reasonable level. Finally, section 7.0 presents the input 
and output for two test cases. 
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4.0 INPUT OPTIONS 


There are two fundamentally different modes of input for the majority 
of the data required by this program. The Doug’ is input mode is very similar 
to the original format described in reference 1. Its main feature is that 
the geometry data is input as a set of points, one or two per card. It is 
therefore fully general, allowing configurations of any geometry whatsoever 
to be input. It is sometimes more convenient to use the Deriv input mode. 

Its main features are the use of namelist input procedures and the inclusion 
of a geometry package of its own, which allows such simple configurations as 
might be used in oarametric studies to be generated by the program with a 
« niinimum amount of input. Choice of the input mode is indicated by a flag 
on the initial carp input. 

Both input modes provide for the following execution options: 


(a) 

Use of the geometry package. 


(b) 

Both lifting and nonlifting components 

present. 

(c) 

Multiple angles of attack. 


(d) 

Planes of symmetry. 


(e) 

Special formulas for the semi-infi ni te 

last elements on wakes 

(f) 

Extra strips. 


(g) 

Ignored elements. 


(h) 

Off-body points. 


(i) 

Intermediate output- 


(j) 

Full or partial execution of the program. 


In addition, the Douglas input mode provides the following options: 

(a) One or two points per card input. 

(b) Piecewise - constant or linear spanwise representation of the 
vorticity on a lifting component. 
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The geometry package provides the following options: 

(a) Transformation of the input coordinates by translation, scaling, 
and/or rotation about an axis. 

(b) Full or partial execution of the geometry oackage. 

(c) Six different algorithms for distributing on-body points on 
N-1 ines. 

(d) Four different algorithms for distributing wake points on N-lines, 

(3) Four different algorithms for distributing N-lines. 

(f) Two modes of operation for distributing N-lines. 

(g) Two different procedures for the final redistribution of points on 
components which intersect other components, 

(h) Three different procedures for the final redistribution of points 
on components which are intersected by other components. 

(i) Punched output at three stages of the calculations. 

Further information regarding the use of the above options can be found 
in references 1, 2, and 3. 
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5.0 INPUT DATA PREPARATIONS 


5.1 Card Input 

Because of the number of options provided, there are many variables 
which can be input to the program and therefore many.tyoes of cards in the 
list below. Most typical cases require only a few of the options, however, 
so a typical deck set-up would be considerably less complicated than shown 
below. The card numbers listed below refer to the types of cards and not 
the actual card count, since some card types are optional and some types 
require more than one card to be input. Multiple cases may be run by 
stacking decks sequentially. 


General Program Cards : 


(a) 

Card #1 

Title card 

Douglas- 

•Mode Input Cards: 

(a) 

Card #1 

Case control card 

(b) 

Card #2 

Uniform onset flow card 

(c) 

Card #3 

Moment origin card (optional) 

(d) 

Card #4 

Component control card 

(e) 

Card #5 

Ignored elements card (optional) 

(f) 

Card #6 

Coordinate cards 

(g) 

Card #7 

Off -body point cards (optional) 

Deriv-Mode Namelist Input Cards: 

(a) 

General 

problem description namelist NPR^B 


5NPR0B 



NB0D=1 , 

NPAN=2, SREF= ... 


SEND 


(b) 

Nonlifting component namelist NB0OY 


SNB0DY 



IB0DY=1 

9 • • • 


SEND 
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5.1 Card Input (Continued) 


(c) Lifting component namelist NPANL 
$NPANL 

IPANEL=1, ... 

SEND 
$ NPANL 

IPANEL=2, ... 

SEND 


Geometry Package Input Cards 


(a) 

Card 

#1 

Control card 

(b) 

Card 

#2 

Component axis card 

(c) 

Card 

#3 

Point distribution card for on-body points 
on N-lines 

(d) 

Card 

H 

Point distribution card for wake points 
on N-lines 

(e) 

Card 


N-line distribution card 

(f) 

Card 

#6 

Multisegment option card 

(g) 

Card 

#7 

Component transformation card 


(optional ) 
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5.2 Input Data Card Descriptions 


The input format for each of the above card types is given in the 
following tables: 

5.2.1 General Program Input Cards 
This card is always input. 


CARD 

TYPE 


REQUIRED 

OR 

OPTIONAL 


CARD 

COLUMN FORMAT 


VARIABLE DESCRIPTION 


1 


Required 1 

2 


9-80 


II 

n 


18A4 


IDERIV Input mcde flag 

IDERIV=0, Deriv input 
IDERIV^O, Douglas input 

IGE0M Geometry package flag 

IGE2IM=0, Skip geometry 
package 

IGE0M=1, Panel isolated 
components only 
IGE0M=2, Panel isolated 
components and 
calculate inter- 
section curves 
IGE0M=3, Panel isolated 
components, cal- 
culate intersec- 
tion curves, and 
repanel intersect- 
ing and intersect- 
ed components. 

TITLE Description of the case to 
be run 
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5.2.2 



i 


Douglas-Mode Input Cards 

These cards are input only if IDERIVj^O on the initial card. 


REQUIRED 

CARD OR CARD 


TYPE OPTIONAL 

COLUMN 

FORMAT 

VARIABLE 

DESCRIPTION 

1 

i 1 Required 

i ' 

1-4 

A4 

CASE 

Case ID (any 4 alphanumer- 
ic characters) 

i 

f 

5-7 

13 

T0TSEC 

Total number of components 
input 

8-10 

13 

M0MENT 

Moment origin flag 
M0MENT=O, Moments calcula- 
ted about the 
origin 

MjfMENT^O, Moment center 
will be input 


11-13 

13 

NPTS 

Coordinate input flag 
NPTS=1 , 1 point per card 
NPTS^l, 2 points per card 

t 

1 

i 

♦ 

14-16 

13 

N0FF 

Off-body points input flag 
Nj!FF=0, no off-body points 
i nput 

N0FFj<O, off-body points 
input 


17-19 

13 

LIST 

Execution flag 
LIST=0, full execution 
LISTjf^O, partial execution 
(basic elements 
formed only) 


20-22 

13 

MPR 

Matrix print flag 

MPR=0, no matrix print 
MPR=1, Vij matrix orint 
MPR=2, Aij matrix & o- solu- 
tion print 

MPR=3, Onset flows &o- solu- 
tion print 

1 

23-25 

13 

I0UT 

Geometry print flag 
I0UT=O, basic geometry print 
I0UT^O, additional geometry 
print 


26-28 

13 

IG 

Ignored elements flag 
IG=0, no ignored elements 
IG^O, ignored elements 
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5.2.2 

Douqlas-Mode Input Cards (Continued) 


CARD 

TYPE 

REQUIRED 

OR 

OPTIONAL 

CARD 

COLUMN 

FORf^T 

VARIABLE 

DESCRIPTION 

1 

Requi red 

29-31 

13 

LASWAK 

Last wake element flag 
LASWAK=0, truncated wake 
LASWAK^O, semi -infinite wake 



32-34 

13 

I AT AC K 

Number of uniform onset flows 
input 



35-37 

13 

I WIDTH 

Spanwise vorticity flag 
IWIDTH=0, piecewise constant 
vorticity 

IWIDTH^O, piecewise linear 
vorticity 



38-40 

13 

IZERO 

Noncirculatory solution flag 
IZER0=O, basic case 
IZER0^O, basic case plus 
noncirculatory 
solution 



41-44 

F4.0 

SYMl 

Symmetry flag for x-z plane 
SYM1»0.0, no symmetry 
SYM1=+1.0, symmetric about 
axis 

SYM1=-1.0, antisymmetric 
about axis 



45-48 

F4.0 

SYM2 

Symmetry flag for x-y plane 
SYM2=0.0, no symmet'^y 
SYM2*+1.0, symmetric about 
axis 

SYM2=-1.0, antisymmetric 


about axis 


2 

Required 1-10 

FIO.O 

ALPHAX 

X-component of the 
uniform onset flow 

first 


11-20 

FIO.O 

ALPHAY 

Y-component of the 
uniform onset flow 

first 


21-30 

FIO.O 

ALPHAZ 

Z-component of the 

first 


uniform onset flow 

For more than 1 uniform on- 
set flow, repeat the same 
format from CC.31-60. Two 
onset flows per card 
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5.2.2 Doug1as-Mode Input Cards (Continued) 


CARD 

TYPE 

REQUIRED 

OR 

OPTIONAL 

CARD 

COLUMN 

FORMAT 

VARIABLE 

DESCRIPTION 

3 

Optional 

1-10 

FIO.O 

JDRI6NX 

X-coordinate of the moment 
center 



11-20 

FIO.O 

pRIGNY 

Y-coordinate of the moment 
center 



21-30 

FIO.O 

pRIGNZ 

Z-coordinate of the moment, 
center 






This card is required if 
MpMENTf^O on card tvpe 1 

4 

Required 

1-4 

14 

NTP(J) 

Lifting component flag 
NTP(J)=0, component J is 
nonlifting 

NTP(J)=1, comoonent J is 
lifting 



5-8 

14 

JSTRIP(J) 

Number of strips in com- 
Donent J 



9-12 

14 

IGG(J) 

Ignored elements flag for 
component J 
IGG(J)=0, no ignored 
elements 

IGG(J)^0, ignored elements 



13-16 

14 

NLINl(J) 

Edge condition flag for 


first strip on component 
J (only needed for lift- 
ing components using the 
piecewise linear vorti- 
city option) 

NLIN1=1, normal strip 
NLIN1=3, continuation 
from previous 
component 

NLIN1=4, extra strip in 
front 
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5.2.2 Douqlas-Mode Input Cards (Continued) 


REQUIRED 

CARD OR CARD 

TYPE OPTIONAL COLUMN FORMAT VARIABLE DESCRIPTION 


4 Required 17-20 14 NLINN(J) 

(Cont) 


21-24 14 IXFLG(J) 


Edge condition flag for 
last strip on component 
J (only needed for lift- 
ing components using the 
piecewise linear vorti- 
city option) 

NLINN=1 , normal strip 
NUNN=2, next comoonent 
is a continua- 
tion of this one 
NLINN=4, extra strip at 
end of component 

Extra strip flag for com- 
ponent J 

IXFlG= 1, first strip is 
an extra strip 
IXFLG=3, last strip is an 
extra strip 

IXFLG=2, first and last 
strips are extra 
strips 

If NTP(v))=0, the last 
three variables on this 
card need not be input 


5 

Optional 1-4 

14 

IG1(I,J) 

First ignored element on 
strip I of component J 


5-8 

14 

IGN(I,J) 

Last ignored element on 


strip I of component J 




6 strips per card, JSTRIP(J) 
total 

This card is required only 
when IGj<0 on card type 1 
and IGG(J)j*0 on card type 4 

Card types 4 and 5 are 
input for all components 
before inout card type 
6 
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5,2.2 


CARD 

TYPE 

6 


Douglas-Mode Input Cards (Continued) 


REQUIRED 

OR 

OPTIONAL 


CARD 

COLUMN 


FORMAT VARIABLE 


DESCRIPTION 


Required 1-10 

11-20 
21-30 

FIO.O 

FIO.O 

FIO.O 

*] 

1 x,y,z coordinates of a 

[ point (all points input 

) N-line by N-line, component 

by component) 

31 

11 

STATUS 

STATUS=0 same N-line 
STATUS=1 new N-line 
STATUS=2 new component 
STATUSES last point input 
(exception below) 
STATUS=4 first point on 
wake 

STATUSES first point on 
wake and also 
last point input 

33-42 

43-52 

53-62 

FIO.O 

FIO.O 

FIO.O 


1 x,y,z coordinates of the 

I next point input 

63 

11 

STATUS 

Same meaning as above 




Repeat cards to input more 
points. If NPTS=1 on card 
type 1 , input only 1 point 
per card (columns 1-32). 




Card type 6 is input for 
all components before input 
of card type 7 


Optional 1-10 

11-21 
21-30 

FIO.O 

FIO.O 

FIO.O 

X0FF 

Y0FF 

Z0FF 

1 x,y,z coordinates of an 

i off -body point 

31 

11 

STATUS STATUS=0, more points input 


STATUS=3, last point input 

32-41 

FIO.O 

X0FF) 

1 

, Same meaning as above 

42-51 

FIO.O 

Y0FFJ 

52-61 

FIO.O 

ZPFF 

1 

62 

11 

stat) 

1 


ORIGINAL PAGE IS 
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5.2.3 Deriv-Mode Input Cards 

These cards are input only if IDERIV=0 on the initial card. 

NAMELIST/NPROB/ - GENERAL PROBLEM DESCRIPTION 
(1) Configuration Data 

NBOD Total number of bodies used to represent the configuration. 

NPAN Total number of panels used to represent the configuration. 

SREF Total configuration reference area. 

CREF Total configuration reference length. 

BREF Total configuration reference span. 

0RIGNX ) 

^RIGNY ? XYZ coordinates of the moment origin 
^RIGNZ ) 

XMACH Mach number (default value 0.). 


( 2 ) 

lATACK 

ALPHA 

BETA 

DEL0SX ] 
DEL0SY > 
DEL0SZ ) 


Onset Flow Data 

Number of angles of attack and/or sideslip (Max 10). 

Array of angles of attack in degrees (default values 0.). 

Array of angles of sideslip in degrees (default values 0.). 

XYZ components of the delta onset flow corresponding to the panel 
centroids (default values 0.). 


(3) Execution Options 


SYMl X-Z plane symmetry flag 

0, no symmetry 

1, symmetrical solution 

-1, anti symmetrical solution 

SYM2 X-Y plane symmetry flag (SYM2/0 only if SYMl/0). 

0, no symmetry 

1, symmetrical solution 

-1, anti symmetrical solution 

LIST Execution flag 

0, full execution 

1, basic elements are formed only to check geometry 

I|8UT If nonzero, then extended geometric quantities printed. 
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5.2.3 Denv-Mode Input Cards (Continued) 


MPR 


N^FF 

X0FF) 

YPFFJ 

Z0FF) 

LASWAK 


IG 


IZER0 


( 1 ) 

IB0DY 

BAREF 

BCREF 

BCTI 


XB0) 

YB(J> 

ZBjj) 


Matrix print flag. 

1, Vij matrix printed 

2, A^j matrix and sigma solution printed 

3, onset flows and sigma solution printed 

Number of off -body points at which velocities calculatad (Max 200). 


XYZ coordinate arrays of the off-body points 


Last wake element flag 
0, truncated wake 
j*0, semi infinite wake 

Ignored elements flag 
0, no ignored elements 
^ 0 , ignored elements 

Noncirculatory solution flag 
0, basic case 

j^O, basic case plus noncirculatory solution 
NAMELIST/NBODY/ - BODY DESCRIPTIVE AND MODELING DATA 


Coordinate System Definition and Symmetry Options 
Number of this body. 

Body reference area. 

Body reference length. 

Coordinate system flag (default value 1). 

0, input in cylindrical coordinates 

1, input in cartesian coordinates 


XYZ coordinates of body coordinate system origins 


(2) Cross-Section Input Data 

NXS Number of stations for cross-section input (Max 41). If zero, the 

body is constructed by transforming a unit radius cylinder according 
to the multiplication and displacement tables below and XS, NTY, TYS, 
and RZS need not be input. 

XS Array of X-stations at which cross-sections input. 

NTY Number of circumferential input points at each X-station (Max 21), 
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5.2. 

TYS 

RZS 


NXM 

XMC 

YM 

ZM 

DY 

DZ 

NVX 

BLjSI 

VX 

NVY 

BTHI 

VY 

IGll 

IGNN 


Deriv-Mode Input Cards (Continued) 

Array of (NTY.NXS) circumferential 0 or Y coordinates. Input is 
clockwise from the top of the body. The e or Y coordinates for 
the first section are input as TYS(1,1), for the second as TYS(1,2), 
etc. 

Array of (NTY.NXS) circumferential R or Z coordinates correspond- 
ing to the TYS and XS stations. The R or Z coordinates for the 
first section are input, then the second, etc. 

(3) Multiplication and Displacement Tables 

Number of X-stations for multiplication and displacement grid (Max 49). 

Array of X-stations for multiplication and displaceinent grid. Omit 
if same as XS array. 

Y multiplication factor array. (Default values 0.) 

Z multiplication factor array. Omit if all cross-sections circular. 
Array of Y displacements 
Array of Z displacements 

(4) Subpanel Modeling Data 

Number of longitudinal subpanels down the body (Max 300). 

Longitudinal subpanel distribution flag (default va''ue 1). 

-1, spacing input in VX array 

0, cosine spacing VXj=BCREF/2[l-cos(j-l )tt/NVX] j-1 ,NVX+1 

1 , evenspacing VXj=BCREF/NVX[j-l] j=l ,NVX+1 

Array of X locations defining longitudinal subpanel edges 
(NVX+1 values input - necessary only if BL)5I=-1). 

Number of circumferential subpanels around the cross-sections (Max 50). 

Circumferential subpanel distribution flag (default value 0). 

0, equal increments of theta around cross-section 

1, theta angle to each subpanel edge, measured in degrees clock- 
wise from the top of the body, input in VY array. 

-1, TYS and RZS geometry arrays define subpanel edges. 

Array of angles defining circumferential subpanel edges 
(NVY+1 values - input only necessary if BTHI=1 ) 

Array of indices of first ignored element on each strip of the body. 

Array of indices of last ignored element on each strip of the body. 
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5.2.3 Deriv-Mode Input Cards (Continued) 


NAMELIST/NPANL./ - PANEL DESCRIPTIVE AND MODELING DATA 

(1) Coordinate System Definition and Symmetry Options 
IPAN Number of this panel 


XP0) 

YPO > XYZ CO. rdinates of panel coordinate system origins. 

m) 


R0TA ) 
R0TB } 
R0TC ) 


First, second, and third rotations in degrees about axes parallel to 
the X- , Y-, and Z-axes, respectively. 


R0TX 1 

R0TY ? XYZ coordinates of the point about which the panel is rotated. 


(2) Planfc mi Description 

NPL tlumber of points describing planform leading edge (Max 20). 

PXL) 

PYL > Array of XYZ coordinates of leading edge points. 

PZL ) 

CH0RD Amy of chord distances at leading edge points. 


(3) Camber and Twist Data 

NET Num er of spanwise stations where camber and/or twist is described. 

If 10 bpanwise variation, set NET=1 (Max 30). 

PET Array of NET spanwise stations (s/as where s»computed semi span). 

NXC Number of chordwise stations where the mean camber line is described 

(Max 30). 

PXC Array of NXC chordwise (x/c) stations. 

CAMBER Array of (NCX.NET) camber values input as fraction chord (z/c). 

Camber values for u'C first spanwise station input as CAMBER(1,1), 
for the second >'■. CAMBER(1,2), etc. 

TWIST Array of NET twist values at the spanwise stations defined in the 
PET arra>.'. (Input in degrees, positive leading edge up.) 

TWPCR Array of NET chordwise (x/c) locations about which airfoil sections 
twi ed. 
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5.2.3 Deriv-Mode Input Cards (Continued) 

(4) Thickness Data 

NETA Number of spanwise stations where thickness input. If no spanwise 
variation, set NETA=1 (Max 30). 

ETAT Array of NETA spanwise stations. 

NXCT Number of chordwise stations at which thickness input (Max 30). 

XjJCT Array of NXCT chordwise (x/c) stations. 

THICK Array of (NXCT, NETA) half- thickness (t/2c) values input in the same 
manner as the camber tables. If the same stations are used for the 
thickness and camber input, NETA, ETAT, NXCT, and X£)CT need not be 
input. 

(5) Subpanel Modeling Data 

NVS Number of spanwise subpanels (Max 40). 

PTHI Subpanel spanwise distribution flag (default value 0). 

0, spacing at even span increments VSj=(j-l )/NVS j=l ,NVS+1 

1, spacing input as fraction semispan in VS array 

-1 , cosine spacing VSj=l-cos[(y-l )(ir/2)/NVS] I 

VS Array of spanwise (s/as) subpanel edges (NVS+1 values - input 

necessary only if PTHI=1). 

NVC Number of chordwise subpanels (Max 40) 

PLOI Subpanel chordwise distribution flag (default value 0). 

-1, spacing at even chord increments VCj=(j-l )/NVC j=l ,NVC+1 

0, cosine spacing VCj=0.5[l-cos(j-l )tt/NVC] 

1, spacing input as fraction unit chord in VC array 

2, spacing with cos/sinh distribution 

VC Array of chordwise (x/c) stations defining suboanel edges 

(NVC+1 values - necessary only if PL0I=1). 

NPWV Number of wake panels (Max 15). 

lEXTRA Extra strip flag. 

1, outboard strip is extra 

2, outboard and inboard strips are extra 

3, inboard strip is extra 

IGll Array of indices of first ignored element on each strip of the panel. 

IGNN Array of indices of last ignored element on each strip of the panel. 
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5.2.4 Geometry Package Input Cards 


on the initial card. 


DESCRIPTION 

Specified number of on-body points 
per N-line on component J. 

Specified number of wake points per 
N-line on component J 

Specified number of N-lines on 
component J 

Spacing algorithm for on-body 
points on N-lines 
NALGB=0, input distribution, 
unaltered 

NALGB=1 , input distribution, 
auqmented 

NALGB=2, constant increments in - 

arc length 

NALGB=3, cosine distribution 

NALGB*4, curvature-dependent 

distribution * 

NALGB=5, user-specified |' 

distribution , 

Spacing algorithm for wake points 

on N-lines. (Values have same 

significance as values of NALGB, 

except that NALGW=3 and NALGW=4 

should not be used.) j 

Spacing algorithm for N-lines ■ r 

NALGS=0, input distribution, ' 

unaltered 

NALGS=1 , input distribution, 

augmented i, 

NALGS=2, constant increments 
NALGS=3, user-specified distri- 
bution i’ 




5.2.4 Geometry Package Input Cards (Continued) 

REQUIRED 

CARD OR CARD 

TYPE OPTIONAL COLUMN FORMAT VARIABLE DESCRIPTION 

1 Required 25-28 14 IM0DE(J) Planar-section/arc-length mode 

(Cont) flag. 

IM0DE=O, planar-section mode 

using default values of 
direction cosines of 
axis [(1,0,0) for non- 
lifting components; 
(0,-1 ,0) for lifting 
components)]. 

IM0DE=1 , planar-section mode 
using input values of 
direction cosines of 
axis. 

IM0DE=2, arc-length mode. 


29-32 

14 

NSEG(J) 

Number of segments into which 
N-lines are borken (default*!). 

33-36 

14 

NTR(J) 

Comoonent transformation flag. 

(The value indicates the number 
of transformations to be performed.) 

37-40 

14 

INTS(J) 

Index of the component (if any) 
which component 0 intersects. 

41-44 

14 

INTD(J) 

Index of the component (if any) 
which intersects component J. 

45-48 

14 

NALG2(J) 

Repaneling flag for intersecting 
components 


NALG2=0, minimal repaneling 
NALG2=1 , full repaneling 


ORIGINAL PAGE IS 
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5.2.4 Geometry Package Input Cards (Continued) 

REQUIRED 

CARD OR CARD 

TYPE OPTIONAL COLUMN FORMAT VARIABLE DESCRIPTION 

1 Reauired 49.52 14 NAL63(J) Repaneling flag for intersected 

(Cort) components 

NALG3*0, no repaneling 
NALG3=1 , (a) nonlifting compo- 
nents. Full repaneling 
with an N-line passing 
through every point on 
intersection curve. 

(b) lifting components. 
Full repanel inq, but 
with gaps in the region 
of the intersection 
curve. 

NALG3=2, (a) nonlifting compo- 
nents. Full repaneling 
with an N-line pasnng 
through every second 
point on the intersec- 
tion curve. 

(b) lifting components. 
Full repai.cling with 
additional nonlifting 
elements to fill the 
gaps around the inter- 
section curve. 

53-56 14 IPl(J) Punch flag for transformed input 

coordinates 
IP1=0, no punch 
IP1?*0, punch 

57-60 14 IP2(J) Punch flag for coordinates after 

repaneling isolated components 
IP2=0, no punch 
IP2^0, punch 

61-64 14 IP3(J) Punch flag for final redistributed 

coordinates 
IP3*0, no punch 
1P3?*0, punch 
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5.2.4 Geometry Package Input Cards (Continued) 


REQUIRED 


CARD 

OR 

CARD 




TYPE 

OPTIONAL 

COLUMN 

FORMAT 

VARIABLE 

DESCRIPTION 

2 

Optional 

1-10 

FIO.O 

C0SX(J) 

X-direction cosine of the axis of 
component J. 



11-20 

FIO.O 

C0SY(J) 

Y-direction cosine of the axis of 
component J. 



21-30 

FIO.O 

C0SZ(J) 

Z-direction cosine of the axis of 
component J. 

This card is required if IM0DE(J)=1 
on card type 1 

3 

Optional 

1-10 

8F10.0 

SSB(J.l) 

Specified arc lengths of on -body 



11-20 


SSB(J,2) 

points on N-lines on component J 



21-30 



0.0 < SSB < 1 " 



31-40 



8 poTnts per card, NB(J) total 



41-50 

51-60 

61-70 

71-80 



This card is required if NALGB(0)=5 
on card type 1 

4 

Optional 

1-10 

8F10.0 

SSW(J.l) 

Specified arc lengths of wake points 


11-20 


SSW(J,2) 

on N-lines of component J 



21-30 



0.0 < SSW < 1.0 



31-40 


• 

8 points per card, NW(J) total 



41-50 

51-60 

61-70 

71-80 


« 

This card is required if NALGW(J)=5 
on card type 1 

5 

Optional 

1-10 

8F10.0 

SSS(J,1) 

Specified distribution of N-lines 


11-20 


SSS(J,2) 

on component J 



21-30 


• 

0.0 < SSS < 1 .0 



31-40 


• 

8 poTnts per card, NS(J) total 



41-50 

51-60 

61-70 

71-80 


• 

This card is required if NALGS(J)=3 
on card type 1 
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5.2.4 Geometry Package Input Cards (Continued) 


REQUIRED 

CARD OR CARD 

TYPE OPTIONAL COLUMN FORMAT VARIABLE DESCRIPTION 


6 Optional 1-4 514 

5-8 
9-12 
13-16 
17-20 


IEND(J,1) Point numbers of the ends of the 
IEN0(j,2) segments on M-lines on component J 
. NSEG(J) total (max = 5) 

. This card is required if NSE6(J)>1 
on card type 1 


7 Optional 1-4 14 ITR2(J,K) 


11-20 FIO.O TR(J,K,1) 


21-30 FIO.O TR(J,K,2) 


31-40 FIO.O TR(J,K,3) 


41-50 FIO.O TR(J,K,4) 


Type of transformation to be 
applied to component J 
ITR2=1, scaling 
ITR2=2, translation 
ITR3=3, rotation 

First transformation parameter: 

X multiplication factor if ITR2=1 
X translation if ITR2=2 

angle of rotation if ITR2«3 

(positive counterclockwise, looking 
in the direction of the axis). 

Second transformation parameter: 
y multiplication factor if ITR2»1 
y translation if ITR2»2 

X direction .osine of axis 
of rotation if ITR2»3 

Third transformation parameter: 
z multiplication factor if ITR2=1 
z translation if ITR2»2 

y direction cosine of axis 
of rotation if ITR2=3 

Fourth transformation parameter: 
Dummy if ITR2*1 or ITR2=2 
z direction cosine of axis of 
rotation if ITR2=3 


Note : The axis of rotation is 
assumed to pass through the origin. 
This card is required if ITR(J)^0 
on card type 1. Up to 10 trans- 
formations may be performed, in any 
order. 
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6.0 OVERLAY STRUCTURE 


In order to minimize the core requirements, this program uses several 
overlays. The general overlay structure and the functions of each overlay 
level are described below. 

. Overlay (0,0) - Program NUED 

. Function - Controls overall program execution 

Subroutines - HEADER, READl , WRITEl , READS, WRITES, BL0CK DATA, 
READRS, WRITRS, RITER, BUFFIN, REWIND, SKIPX, ZER0, 
C0DMT, C0DIM, MI SI 

Overlay (1,0) - Program INITIAL (IDERIV=1) 

Function - Data input using Douglas mode 
Subroutines - INPUT 

. Overlay (10,0) - Program READ (IDERIV=0) 

. Function - Data input using Deriv mode 

Subroutines - SQUEEZE, B0DYIN, B0DGE0M, PANELIN, B0DSYM, SM00TH, 
CHLSKY, PANPRT, GPRINT, B0DPRT 

Overlay (12,0) - Program V0RB0D (IDERIV=0) 

. Function - Generates body geometry data 

Subroutines - XYZVB, MQPTS, TRPS 

. Overlay (1S,0) - Program V0RPAN (IDERIV=0) 

Function - Generates wing geometry data 

Subroutines - XYZVW, PANDAT, LINIM, PERIM, PSL0PE, LFTDRV, 

LINE, VLINE 

Overlay (7,0) - Program REPAN (IGE0M>O) 

. Function - Controls geometry package calculations 

Subroutines - INTRP2, DERV, INTERP, ARCLGN, DERV2, CUBPLN, 

N0RM, AMAXA 

. Overlay (7,1) - Program GEINPT (IGE0M>O) 

. Function - Data input for geometry package 
, Subroutines - TRANSF 

MCOOMMnLM. DOUaUAS ,onf»onATioM 
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Overlay (7,2) - Program PANELC (I6E0M^O) 

. Function - Augmentation and redistribution of points on N-lines 

. Subroutines - SPACEC 

Overlay (7,3) - Program PANELS (IGE0M>O) 

. Function - Augmentation and redistribution of N-lines 
Subroutines - SPACES, CBPLN2 
Overlay (7,4) - Program CR0SS (IGE0M^2) 

. Function - Calculates geometric cross-derivative terms 

. Subroutines - None 

Overlay (7,5) - Program INTSEC (IGE0M^2) 

. Function - Calculates intersection curves between comoonents 

Subroutines - AMINA,AMAX3, PLANE, B0UNDS, SIDES, REPLN, UWPLN, 
PC0EFS, CUBCUB 

Overlay (7,6) - Program REPANl (IGE0M>2) 

. Function - Repanels intersecting components 

. Subroutines - INTRP3 

Overlay (7,7) - Program REPAN2 (IGE0M>2) 

. Function - Repanels intersected components 

Subroutines - REPAN3, PRN3C, SHUFL, REPAN4, SRCHl 
Overlay (20,0) - Program ELMNTS 

Function - Controls calculation of element geometric quantities 
Subroutines - DNLIFT, DLIFT 
Overlay (2,0) - P; ^gram VF0RM 

Function - Controls velocity matrix formation 
. Subroutines - None 

Overlay (2,1) - Program VFMNLF 

. Function - Forms velocity matrices for nonlifting components 
Subroutines - None 
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Overlay (2,2) - Program VFMLFT 

. Function - Forms velocity matrices for liftinq components 

. Subroutines - NEAR, WNEAR 

Overlay (3,0) - Program PSWISE (IWIDTH ^ 0 ) 

. Function - Computed onset flows for piecewise linear spanwise 
vorticity option 

. Subroutines - None 

Overlay (4,0) - Program AF0RM 

. Function - Computes normal velocity matrix 

Subroitines - None 

Overlay (5,0) - Program C0LS0L 

. Function - Solves the system of linear equations 
. Subroutines - REQFL 

Overlay (6,0) - Program C0MFL0 

. Function - Controls final calculations and output 

. Subroutines - None 

Overlay (6,1) - Program PKUTTA 
. Function - Combines individual flow solutions 

Subroutines - None 
Overlay (6,2) - Program VELPRS 

Function - Calculates final velocities 
. Subroutines - PRINT 
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7.0 TEST CASES 


The following section describes two test cases which can be used to 
verity the execution of the program. 

7.1 Test Case Number 1 

The first test case included in this report is a simple swept, tapered 
wing input using the Deriv-mode option. The planform is specified by two 
points on the leading edge (at the root and tip) and two chord values. The 
cross-section shapes are specified by a table of the thickness distribution. 
There is no camber or twif*". The geometry package is used (in its simplest 
mode) to produce a wing paneling consisting of three spanwise strips (four 
N-lines), ten on-body elements (eleven points) per strip, and one wake element 
per strip. The potential -flow pressure distribution is calculated using the 
options for piecewise-constant spanwise vorticity distributions, semi -infinite 
last wake elements, and one plane of symmetry. 
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7.1.1 


Input Data 


Card 

Column 

I ? 3 4 « 6 7 

12 34 56 7BOO I ? 345*- 7B9012 34567P9012 34 5678901234567P*J01 234 567 890 12 34 5678 90 

I t^FOT TflPFRFO WIMG (n«=RIV IMPUT) 

iNPQPB 

NPaN=i,5PFP = l.3.C»EF=l. Jf 9RFF = l.D, I4TarK=2,flLPMa=13. ,3., 

*:vmi = 1.0,NOFF=2,XPFF=-.U-.5, Yr>FF = .5 , . 5 ♦ 7PF F = .0, .0.LASWAK= I , 

3 FNP I 


JNPavL 

IPANFL=l.N''»WV*l ,N*PL = 2i'''eTA=l,NXCT = 6fNVS = l,Nvr = 5, 
PXt. =0., l.760e2«PYL=.0,4.,P7l = . J,. 3» 

CHpffp=l .3P2319, O. 69066 » 

XOr T=0. ..006294,.031401,.351299,.746699,1.0, 
thick =0., .006 730» .014060, .0375 70, . C22 1 3 0, . 00 , 
SFNO 


11 1430 ? 
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7.1.2 Output Data 

The output for the first test case is given in the following pages. 
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7.2 Test Case Number 2 


The second test case Is a swept, tapered wing intersecting a fuselage con- 
sisting of a cylindrical midsection with ellipsoidal fairings fore and aft. The 
Douglas-mode input option is used. The shape of the wing is specified by input- 
ting points on the tip and root N-lines. The shape of the fuselage is specified 
by inputting points on ten N-lines. Since the midsection of the fuselage is 
cylindrical, the fuselage is paneled using the multisection component option of 
the geometry package. Both wing and fuselage are paneled using the planar- 
section mode of distributing N-lines. The intersection curve is calculated and 
the wing and fuselage are repaneled. The final element distribution nn the wing 
consists of three strips (including one extra strip at the root), four on-body 
elements per strip, and one wake element per strip. The final element distribu- 
tion on the fuselage consists of two strips forward of the wing, two strips aft 
of the wing, one strip above the wing, and one strip below the wing. There are 
four elements per strip on the portions forward and aft of the wing and two 
elements per strip on the portions ab «ve and below the wing. The potential- flow 
pressure distribution is calculated using the options for piecewise-constant 
spanwise vorticity distributions, semi -infinite last wake elements, one plane 
of symmetry, and one extra strip on the wing. 
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7.2.1 


Input Data 


Card 

Column 

12 3 4 5 6 

1234567890 I 234567890 123456 T 890123456780012345678901234567890 

13 WlMr.-FUSELAGE TEST CASE 



10 

30 


1 0 
0.0 0.0 
3 


1 


2.451475 

4.000271 

0.0 2 

2.276530 

4.000271 

- 0.0152840 

2.003443 

4.000271 

- 0.0259480 

1.782502 

4.000271 

- 0.0097110 

1.765165 

4.000271 

- 0.0046480 

1.760815 

4.000271 

0.0 0 

1.765165 

4.000271 

0.0046480 

1.742502 

4 . U 0 U 271 

0.0097110 

2.003443 

4.000271 

0.0259480 

2 , Z 1 b 53 Q 

4.000271 

0.0152840 

2.451475 

4.000271 

0.0 0 

2 . 5377^0 

4.000000 

0.0 4 

1.382319 

0.0 

0.0 1 

1.032176 

O.O 

- 0.0305910 

0.485608 

0.0 

- 0.0519340 

0.043406 

0.0 

- 0.0194350 

0.008710 

0.0 

- 0.0093030 

0.0 

0.0 

U.O 0 

0.008710 

0.0 

0.0093030 

0.043406 

0.0 

0.0194350 

0.435608 

0.0 

0.0519340 

1.032176 

0.0 

0.0305910 

1.382319 

0.0 

0.0 0 

1.554960 

0.0 

0.0 4 

- 3.333333 

0.0 

0.0 2 

- 3.333333 

0.0 

0.0 0 

- 3.333333 

0.0 

0.0 0 

- 3.333333 

0.0 

0.0 0 

- 3.333333 

u.o 

0.0 0 

- 3.333333 

0.0 

0.0 0 

- 3.333333 

0.0 

0.0 0 

- 3.250000 

0.0 

0.1425461 

- 3.250000 

0.071273 

0.1234490 

- 3.250000 

0.109197 

0.0916270 

- 3.250000 

0. 142546 

0.0 0 

- 3.250000 

0.109197 

- 0.0916270 

- 3.250000 

0.071273 

- 0.1234490 

- 3.250000 

0.0 

- 0.1425460 

- 3.000000 

0.0 

0.2762601 

- 3.000000 

0.138130 

0.2392480 

- 3.000000 

0.211627 

0 . 177576 C 

- 3.000000 

0.276260 

0.0 0 

- 3.000000 

0.211627 

- 0 . 177576 C 

- 3.000000 

0.138130 

- 0.2392480 

- 3.000000 

0.0 

- 0.2762600 


1.0 
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7.2.1 (Continued) 

Card 

Coliflnn 

12 3^56 

123456780012345678901234567890123456789012345678901234567800 



-2.600000 

0.0 

0.3866711 


-2.600000 

0.193436 

0.335040G 


-2.600000 

0.296361 

0.2486760 

■ - 

-2.600000 

0.386871 

0.0 Q 


-2.600000 

0.296361 

-O.240606C 


-2.600000 

0.193436 

-0.3350400 

s 

-2.600000 

0.0 

-0.3868710 


-1.333333 

0.0 

0.5000001 

‘I 

' -1.33333? 

U.250OU0 

0.4330130 


-1.333333 

0.383022 

0,3213940 

J 

1 -1.333333 

0.500000 

0.0 0 


-1.333333 

0.363022 

-0.3213940 


-1.333333 

0. 250000 

-0.A330130 


-1.333333 

0.0 

-0.5COOOOO 

^ 2.400000 

0.0 

0.5000001 

2.400000 

0.250000 

0.4330130 


2.400000 

0.383022 

0.3213940 


2.4U0C00 

0.5UOOOU 

0.0 0 


2.400000 

0.383022 

-0.3213940 


2.400000 

0.250000 

-0.4330130 

Sr-4|P 

2.400000 

0.0 

>0.5000000 

•<** 

•r>r« 

3.500000 

0.0 

0.4545301 


3.500000 

0.227265 

0.393634C 


3.500O0O 

0.348190 

0. 2921660 


3.500000 

0.454530 

0.0 0 


3.500000 

0.348190 

-0,2921660 


3.500000 

0.227265 

-0.3936340 


3.500000 

0.0 

-0.454530C 


4.300000 

0.0 

0.2035551 

$ 

4.30000U 

0.144277 

U. 2498960 

?rV 

4.300000 

0.221046 

0.1854600 


4.300000 

0.288555 

0.0 0 


4.300000 

0.221046 

-0. 1854800 


4.300000 

0.144277 

-0.249896C 

• 

■j6^ 

4,300000 

0.0 

-0.28R5550 

4 

4.500000 

O.O 

0.1998261 


4.500000 

0.099913 

0.1730550 

A«i 4.500000 

0.153076 

0. 1284460 

«* 4.500000 

0.199826 

0.0 C 

4.500000 

0.153076 

-0.1284460 


4.50000) 

0.099913 

-0.1730550 


4.500000 

g 

0.0 

-0. 1998260 

s 

1 4.666667 

0.0 

0.0 1 


1 4.666667 

0.0 

0.0 0 


■ 4.666667 

0.0 

0.0 0 


4.666667 

0.0 

0.0 0 

_ 

} 4,666667 

0.0 

0.0 0 


1 4.666667 

0.0 

0.0 0 


4.666667 

j 

0.0 

O.O 3 


I 5 1 

3 3 0 

2 


5 0 

5 2 0 

2 


5 6 

10 



2 
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7,2.2 Output Data 


The output for the second test case Is given In «.. < following pages. 
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